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Abstract 

The Einstein-Klein-Gordon field equations arc solved in a inhomo- 
geneous shear-free universe containing a material fluid, a self-interacting 
scalar field, a variable cosmological term, and a heat flux. A quintessence- 
dominated scenario arises with a power-law accelerated expansion com- 
patible with the currently observed homogeneous universe. 



1 Introduction 

Recently, there have been claims in the literature that the Universe, besides 
its content in normal matter and radiation, must possess a not yet identified 
component (usually called quintessence matter, Q-matter for short) [P, 

characterized by a negative pressure, and possibly a cosmological term. 
These claims were prompted at the realization that the clustered matter 
component can account at most for one third of the critical density. There- 
fore, an additional "soft" (i.e. non-clustered) component is needed if the 
critical density predicted by many inflationary models is to be achieved. 

Very often the geometry of the proposed models is very simple, just 
Friedmann-Lemaitre-Robertson- Walker (FLRW) . In constrast to FLRW mod- 
els, inhomogeneous spaces are in general compatible with heat fluxes, and 
these might imply important consequences such as inflation or the avoidance 
of the initial singularity Q. Here we focus on an isotropic but inhomoge- 
neous spherically symmetric universe which besides a material fluid contains 
a self-interacting scalar field (which can be interpreted as Q-matter), and a 
cosmological term, A which, in general, may vary with time. 
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Density inhomogenities triggered by gravitational instability must be 
present at any stage of evolution. We only mention that the negative pres- 
sure associated to Q-matter and A will tend to slow down the growing modes 
(see e.g. P), and shift the epoch of matter-radiation equality toward 
more recent times. 



2 Einstein-Klein-Gordon field equations 

Let us consider a shear-free spherically-symmetric spacetime with metric 

i 

-v{t,rf df + dr^ + r^dl^^j. (1) 

where as usual dO^ = dO"^ + sin^ 0d(j)'^. Units have been chosen so that 
c = G = 1. As sources of the gravitational field we take: a fluid of material 
energy density pj = pf{r,t), hydrostatic pressure pj = Pf{r,t), with a radial 
heat flow (q^ = qrir,t) and qt = Qe = Qcp = 0), plus a cosmological term, 
related to the energy density of vacuum by A = Sirp^ac, that depends only on 
time A = A{t), and a self-interacting scalar field (f) driven by the potential 
V{(p) whose equation of state is = (7,^ — 1) p,^. Hence the scalar field 
can be interpreted as Q-matter -see e.g. |^]. The stress energy-tensor of 
the normal matter, with a heat flow, plus Q-matter (scalar field) and the 
cosmological term is 




U = (Pf +P<l>+Pf+ P<t>)u'uk + {A-pf- p^)6i + q'uk + qku' , (2) 

As equation of state for the fluid we choose pf = {jf — 1) pf where 7/ is a 
function of t and r. Taking into account the additivity of the stress-energy 
tensor it makes sense to consider an effective perfect fluid description with 
equation of state p = {'j — 1) p where p = Pf + p^j,, p = pf + p^ and 

^ = JiPl±JiPl, (3) 

is the overall (i.e. effective) adiabatic index. The requirement that the 
cosmological term A is just a function of t leads to the restriction that 7 
also depends only on t to render the system of Einstein equations integrable. 
The nice result we are seeking is a solution that has an asymptotic FLRW 
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stage, with A evolving towards a constant, and the heat flow vanishing in 
that hmit 0. 

To write the Einstein equations we use the ansatz F = a{t) + b{t)x and 
V = c{t) + d{t)x with the constraints a{t) d{t) = b{t) c{t). This set of metrics 
contains those of Modak (6 = 0) M, Bergmann (c = a,d = b) @ and Maiti 
[b = d = ka/A, with A; = 0, ±1) ||T^. Another possibility arises when d = 0, 
then re-defining the time by vdt dt, the Einstein equations are 

/) + A = 12a6 + 3d2 +6dte + 362a;2, (4) 



p-A = (2bb - 36^) + 2 (26^ - 3db + ab + db^ x - Sab - Sct^ + 2ad , (5) 

Qr = —4:^/xb {a + bx)'^ , (6) 

where x = and the overdot indicates d/dt. Imposing that A = A(i), the 
general solution to these equations has the form 

a = —2 exp dt ivj J dtw'^ J — ^ , b = exp (^J dt vj^ , (7) 

where w = 2/ J dt{2 — 37), provided 7 / 2/3. Inserting (Q) in @j, ^ and 
(^) we easily compute the cosmological constant and the heat flow. Finally 
the FLRW metric is 



ds 
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-dr^ + R^ (dr^ + r^dJ^^^] , (8) 



where we have introduced the time coordinate dr = dt/a, M = b/a and 
R = 1/1 a|. This metric is conformal to FLRW, and the conformal factor 
approaches unity when M ^ 0. 

3 Constant adiabatic index 

When 7 is a constant different from 2/3, the general solution of @ and @ 
becomes 

a{t) = CiAt~^ + C2At~^ - -KAt^^ , (9) 

3 
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M{t) = k{Ci + CsAt-i - -KAi^ 



(10) 



Two alternatives of asymptotically expanding universes appear depending 
on the map between t and r. 

Case At ^ 



In this limit we obtain 



R{r) 



Co 



A(r) 



2C2 (1 - 37) 
2-37 

3(2-3;)' 



At 



37 

2(37-1) 



4(1-37)'At2' 



, , 8KC2 
Qr(,r,T) ~ ^ r 



2C2 (1 - 37) 



At 



97 
2(1-37) 



(11) 



(12) 



(13) 



2-37 

for 7 7^ 1/3. When 1/3 < 7 < 2/3 we have, for large cosmological time 
T, an accelerating universe that homogenizes with vanishing cosmological 
term and heat flow. In this stage we have a final power-law expansion 
era. For 7 = 1/3 we have asymptotically a de Sitter universe with finite 
limit cosmological term. For the remaining values of 7 the universe begins 
at a homogeneous singularity with a divergent cosmological term. When 
7 < 1/3, the heat flux asymptotically vanishes near the singularity, while 
for 7 > 2/3 it diverges. 

Case At — 00 



In this limit we obtain 



R{t) ~ - 



K 



K{4- 37) 
3(2-37) 



At 



6(1-7) 
4-37 



A(r) ~ 



24(2-37)^ 



,(r, r) ~ -24 



(4 - 37)^ Ar2 ' 
(2-37)^ r 



(4 - 37)^ Ar3 ' 



(14) 
(15) 
(16) 
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for 7 7^ 2/3. When 2/3 < 7 < 1 the universe homogenizes for large cosmo- 
logical time with vanishing cosmological term and heat flow. When 7 = 1, 
the late time evolution changes to an asymptotically Minkowski stage. For 
1 < 7 < 4/3 the universe starts homogeneously in the remote past with a 
vanishing scale factor, cosmological term and heat flow. For the remaining 
values of 7 the universe begins at a homogeneous singularity with a divergent 
cosmological term. 

An exact solution with explicit dependence on the asymptotic cosmolog- 
ical time r can be found when the integration constants Ci and C2 vanish. 
In such a case the metric is 



n 2-37 

1 + mAr 4-37 



(17) 



where m is a redefinition of the old integration constant K, the adiabatic 
index 7 and rg. The last constant was introduced by scaling the radial 
coordinate r — > ro r. 



4 Asymptotic evolution to a quintessence— dominated 
era 

As a first stage of towards more general scenarios with a slowly time-varying 
7, we will explore a model that evolves towards an asymptotic FLRW regime 
dominated by Q-matter (i.e. the scalar field). We will show that this system 
approaches to the constant 7 solutions for large times found above. In this 
regime the equations of Einstein-Klein-Gordon become 

3H^ c^Pf + ^^^ + V{^)+A, (18) 

4> + 3H^ + ^^-J^c^O, (19) 
d<p 

where H = R/R and a dot means d/dr in this section. In last section we 
found that the general asymptotic solution for the scale factor R{t) oc At" 
has the power-law behaviors (pAj), ( p!^ ) for any value of the effective adiabatic 
index 7. Then, using these expressions and (H) and (ID together with (|l|) 
and dl^), we can investigate the asymptotic limit in which the energy of 
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the scalar field dominates over the contribution of the perfect fluid. In the 
regime that 3a^f > 2 the adiabatic scalar field index can be approximated 
by 



2 

3a 



1+1 



37/ 



where a = Pf/ptf, <^ 1- Inserting these equations in (| 
correction to the effective adiabatic index 



(20) 



we obtain the first 
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1±a/^(37/ 



2) a 



(21) 



The negative branch of (^) yields a consistent asymptotic solution for the 
range | < 7 < |. We note that this solution describes a deflationary stage 
with a limiting exponent a = 1. 

Oftenly power-law evolution of the scale factor is associated with loga- 
rithmic dependence of the scalar fleld on proper time [^]. Thus, assuming 
that 0(t) ~ Clnr with the constant C to be determined by the system of 
equations (|^) an d ([T9| ) , and using these expressions together with ([l^ ) and 
( |l5|) in ( |l8|) and (|19D it follows that the leading term of V{(j)) for large 4> is 



(22) 
= 2, 



Using the dominant value of the effective adiabatic index we flnd 
Vb = 2 and C = l/\/2. 

The models considered in this section are based on the notion of "late 
time dominating fleld" (LTDF), a form of quintessence in which the field 
(f) rolls down a potential V{(f)) according to an attractor solution to the 
equations of motion. The ratio a of the background fluid to the fleld energy 
changes steadily as (p proceeds down its path. This is desirable because in 
that way the Q-matter ultimately dominates the energy density and drives 
the universe toward an accelerated expansion |]ll|, IT! 



5 Concluding remarks 

We have investigated a class of solutions of the Einstein field equations with 
a variable cosmological term, heat flow and a fluid with variable adiabatic 
index that includes those of Modak, Bergmann and Maiti and contains a new 
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exact conformally flat solution. We have found that asymptotically expand- 
ing universes occur when 1/3 < 7 < 1 that homogenizes for large cosmologi- 
cal time with vanishing cosmological term and heat flow. For 1/3 < 7 < 2/3 
the evolution is given by (|lT| ) and corresponds to a power-law accelerated ex- 
pansion for large cosmological time r. On the other hand, when 2/3 < 7 < 1 
even though an asymptotic negative cosmological term occurs, the universe 
evolves toward a decelerated expansion. The particular case 7 = 1/3 leads 
asymptotically to a de Sitter universe with a finite limit for A. We have 
shown that homogeneization occurs also for a time dependent adiabatic in- 
dex provided it has a constant limit, for i ^ and t 00, repectively, and 
is analytic about these points. 

We have carried out a detailed analysis of a model in which Q-matter 
dominates over cold dark matter. This LTDF solution is an attractor be- 
cause, even for large initial inhomogeneities and a wide range of initial values 
for (p and <j), the evolution approaches a common path. It was shown that 
this model can be realized for a wide range of potentials provided they have 
an exponential tail. Our LTDF model only requires that the potential has an 
asymptotic exponential shape for large (f). So, the interesting and significant 
features of our model are: (a) a wide range of initial conditions are drawn 
towards a common evolution; (b) the LTDF maintain some finite difference 
in the equation of state such that the field energy eventually dominates and 
the universe enters a period of acceleration. 
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